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ABSTRACT  ;  -  . 

< 

Formulae  are  derived  to  calculate  the  mass  median  (D0)  and  mass 
mean  (Dj^)  diameters  of  droplets  in  sprays.  These  formulae  are  based  on  an 
inverse  exponential  relationship  which  is  empirically  observed  to  occur 
in  liquid  sprays.  The  median  diameter  divides  the  total  mass  of  liquid  in 
half  while  the  mean  diameter  represents  a  more  commonly  used  statistical 
average  for  identifying  random  samples.  The  formulae  are  used  to  describe 
the  characteristics  of  a  typical  set  of  experimental  drop  data  and  are 
compared  to  conventional  formulae  for  calculating  these  parameters. 
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INTRODUCTION 

When  describing  the  characteristics  of  a  liquid  spray,  two 
of  the  more  commonly  used  parameters  are  the  mass  median  (Do)  and 
mass  mean  (Dm)  diameters  of  the  droplets.  The  mass  median  diameter  re¬ 
presents  the  droplet  diameter  at  which  half  the  mass  of  liquid  is  above 
and  half  the  mass  of  liquid  is  below.  However,  because  the  mass  is  re¬ 
lated  to  the  cube  of  the  diameter,  it  is  found  that  the  mass  median 
usually  lies  close  to  the  largest  sized  drops  observed  on  a  trial.  The 
mass  mean  diameter,  on  the  other  hand,  represents  the  value  for  which  the 
sum  of  the  deviations  in  mass  diameter  from  it  is  zero.  This  is  a  more 
commonly  used  type  of  averaging  term  and  for  drop  diameters  does  not 
rest  as  close  to  the  largest  drop  as  does  the  mass  median  diameter. 

The  model  used  here  depends  upon  an  exponential  relationship 
between  the  cumulative  number  of  drops  and  the  drop  diameter.  This 
relationship  has  been  empirically  observed  to  exist  for  many  different 
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spraying  systems.  The  model  is  used  to  determine  both  D0  and  for  a 
liquid  spray  generated  by  an  explosive  dissemination  and  the  results  are 
compared  to  the  standard  methods  for  calculating  these  parameters. 

THEORY 

The  standard  method  for  determining  the  mass  mean  diameter  or 
mass  median  diameter  is  simply  taken  to  be  the  conmon  statistical  methods 
for  calculating  means  or  medians  (1,2).  Therefore,  the  following  for¬ 
mulae  are  presented  with  no  derivations: 


D0  =  D  when  m  =  Mj/2.0 

(1) 

and 

3  t£1). j  ffjT"1 

°"=  /  Ny 

(2) 

where  D  =  arbitrary  drop  diameter  (mm) 
m  =  arbitrary  drop  mass  (g) 

Mj  =  total  mass  of  all  sampled  drops  (g) 
nj  =  number  of  drops  of  diameter  D-j 
Ny  =  total  number  of  sampled  drops 

To  derive  the  equations  for  the  modified  method  of  evaluating 
D0  and  D^,  the  drop  numbers  and  drop  mass  in  each  diameter  class  are  summed 
from  the  maximum  diameter  to  the  smaller  diameters  using  the  following 
definitions: 


Nk  = 

T.  n-j 

(3) 

’  =  nmax 

Mk  - 

k 

l  m-j 

(4) 

1  =  %ax 
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where  mj  =  mass  of  drops  at  diameter  (g) 

The  observed  drop  number-diameter  distribution  (N  vs  D)  empirically  al¬ 
ways  appears  to  be  well  represented  by  one  of  the  following  two  equations 
(these  equations  are  similar  to  those  of  Rosin  and  Ramler,  Reference  3): 
In  N  =  a  -  bD  (5) 

In  N  =  a  -  bD2  (6) 

where  a  is  the  intercept  at  the  ordinate 
ea  is  the  total  number  of  drops,  Nj 
b  is  the  absolute  value  of  the  slope  of  the  line. 

For  each  drop  class,  assuming  that  the  drops  are  spherical,  it 
is  clear  that: 


"k  -  r  Dk3  nk  <7> 

However,  Equation  3  and  4  imply  that  m^  =  AM^  and  nk  =  ANk  which  when 
substituted  in  Equation  7  give: 

■  f-  °k3  <« 

Allowing  the  class  interval  to  approach  zero,  we  obtain  the  differential 
equation  corresponding  to  Equation  (8): 

dM  =  ££-  D3dN  (9) 

In  order  to  eliminate  N  in  Equation  9,  Equation  5  is  manipulated  to  give: 

dN  =  -bea"b0dD  (10) 

Substituting  10  in  9  gives: 

dM  =  D3e"bDdD  (11) 


Integration  of  Equation  11  is  performed  by  reversing  the  order  of  inte 
gration  (the  order  of  summation  as  defined  in  Equation  4)  to  give: 


UNCLASSIFIED 


UNCLASSIFIED 


4 


=  IpMI  f°  D3e"bDdD 
°k 

npNl  '  r  (bD)3e”bDd(bD) 


/Dk(bD)3e'bDd(bD) 


02) 

03) 


Mk  =  T 

r 

or:  *  *  e&r  | 

The  first  integral  is  the  gamma  function  of  4,  f4.  Solutions  of  this  in¬ 
tegral  are  listed  in  many  handbooks  (Reference  4  for  instance).  The  total 
mass  is  determined  when  Dk  0  or: 


M  *PNT 

MT  =  6b^ 

nr  ffP^T 

ft  ‘  b3 

(14) 

Also: 

1 

"l  -  I  (bDk,  3)] 

(15) 

where: 

KbDfc,  3) 

/°k(bD)3e'bDd(bD) 

Q  (bD)3e'bDd(bD) 

(16) 

values 

of  which  are 

listed  by  Pearson  (5). 

The  mass  mean  diameter  can  be  alternatively  stated  as: 
D 

“  VTTP  Nj 


(17) 


Substituting  14  into  17,  a  new  expression  for  is  obtained: 

■  '•«»/»  <'8> 

The  mass  median  diameter,  D0  is  obtained  from  Equation  (15)  when: 

VMj  =  0.5  =  1  -  I(bD0,  3)  (19) 

Using  Pearson's  Tables  (5),  the  following  values  are  obtained: 
bD0  =  3.671 

or  =  3 . 671 /b  (20) 

Returning  to  Equation  (10)  and  using  Equation  (6)  instead  of  (5),  gives: 

dN  =  -2bDea"bD2dO  (21) 

Substituting  21  into  9  and  reversing  the  order  of  integration  as  before, 
yields: 
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Mfc  =  *pfb  I!  D4e’b°2dD 

j  uk 

(22) 

or:  Mk  =  np^|  f  (bD2)1^  e'^bD2)- 

K  gb15  \  0 

/^(bD2)1-5  e'bC)2d(bD2) 
J  0 

Proceeding  as  before: 

irpNj  r .  0.221 6-rrp  Nt 

*  *  6b^  i2-5  ■  bis 

(24) 

\  =  Mj  [1  -  I ( bD^ , 1  -  5 ) j 

(25) 

„  ^3^7329"  1.0995 

"  V  b1-5  b0*5 

(26) 

D0  =  I  ( bDg  ,1.5)  =  0.5 

(27) 

which  using  Reference  5  yields: 

D0  =  1.477/Vb"1 

(28) 

Equations  18  and  26,  20  and  28  represent  the  mass  mean  and  mass  median 
diameters  that  were  sought.  It  can  be  seen  that  these  are  determined 
simply  from  the  slope  of  the  line,  b,  in  Equations  5  and  6  representing 
the  number-diameter  distributions.  In  practice,  logarithms  to  the  base 
10  are  often  used  to  evaluate  b.  The  slope  obtained  is  of  course  less 
than  when  natural  logs  are  used  by  a  factor  of  2.3026  and  the  equations 
become : 

Dm  =  -0.7890/b  (29) 

D0  =  -1.594/b  (30) 

for  In  N  =  D 

and  Dm  =  0.7245/ /tT'  (31) 

Dq  =  0.9733/  \ IT1  (32) 

for  In  N  *  D2 

It  is  interesting  to  note  that  in  both  cases,  the  mass  mean  diameter  is 
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simply  a  constant  fraction  times  the  mass  median  diameter  (0.495  and 
0.744).  This  results  from  Equations  5  and  6  where  all  the  drops  are  re¬ 
lated  to  all  the  diameters  by  one  equation. 

EXPERIMENTAL 

To  verify  these  expressions,  experimental  data  were  obtained 
from  a  recent  spray  trial.  In  this  trial,  liquid  droplets  were  generated 
by  an  explosive  dissemination.  The  droplets  were  released  at  9  meters 
elevation  and  were  collected  downwind  on  3-way  detector  paper  (NSN  6665- 
21-845-8613).  The  sizes  of  stains  were  related  to  the  corresponding  drop 
diameters  by  a  previous  calibration  procedure  not  discussed  here.  Table 
1  lists  the  drop  diameters  measured  and  the  number  of  drops  of  each  size. 

RESULTS  AND  DISCUSSION 

The  experimental  data  were  fed  directly  to  the  computer  and 
analyzed  by  the  computer  programs  listed  in  Appendix  A.  Table  1  lists 
the  reduced  data  required  for  evaluating  the  averaged  diameters  as  given 
by  Equations  1,  2,  29,  30,  31  and  32. 

The  mass  median  and  mass  mean  diameters  are  listed  in  Table  2 
for  each  of  the  techniques.  Figures  1  and  2  show  plots  of  the  cumulative 
frequency  versus  the  drop  diameter  and  drop  diameter  squared  respectively. 
It  can  be  seen  from  the  figures  that  the  exponential  relationship  bet¬ 
ween  the  cumulative  number  and  drop  diameter  is  a  much  better  approxi¬ 
mation  to  this  experimental  data  than  cumulative  number  vs  drop  diameter 
squared.  This  is  quantitatively  expressed  by  a  significantly  better  re¬ 
gression  coefficient  of  0.990  vs  0.971  respectively. 
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Comparing  the  averaged  diameter  as  calculated  by  Log  N  vs  D 
and  the  standard  method,  it  may  be  seen  that  the  mass  mean  diameters  are 
very  close  (0.386  vs  0.385  respecti vely ) .  However  the  mass  median  dia¬ 
meters  are  significantly  different  (0.781  vs  0.663  respectively).  This 
is  most  likely  due  to  the  effect  of  large  diameter  drops  on  this  para¬ 
meter.  The  large  diameter  drops  contain  so  much  of  the  liquid  that  only 
about  3%  to  5"  of  all  the  drops  lie  above  this  parameter.  Because  they 
are  so  few  in  number,  large  sampling  errors  are  likely  to  occur  and  the 
standard  method  considers  only  those  that  are  sampled.  By  smoothing  the 
data  with  a  straight  line,  however  this  problem  is  circumvented  and  the 
slope  technique  presented  by  this  work  is  a  much  better  method  of  evalua¬ 
ting  mass  mean  or  mass  median  diameters  than  the  standard  summation  method. 

CONCLUSIONS 

An  improved  method  has  been  explained  for  evaluating  the  mass 
mean  and  mass  median  diameters  of  liquid  sprays.  The  standard  summation 
method  has  been  compared  to  this  technique  using  experimental  data  from 
a  typical  spray  trial.  The  technique  was  at  least  as  simple  as  the 
standard  method  and  was  more  accurate  for  the  data  tested. 
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TABLE  U 

TYPICAL  SPRAY 

TRIAL  DATA 

Drop  Diameter, 

D  (mm) 

No.  of 
Drops  n_ 

Accumulated 
Number,  N 

D-’ 

_(mm )  ?_  _ 

D3 

(mm) 3 

nD3 
(mm) 3 

1.165 

1 

1 

1  .357 

1 .581 

0.104 

0.993 

2 

3 

0.987 

0.981 

1.961 

0.917 

1 

4 

0.841 

0.772 

0.772 

0.898 

1 

5 

0.807 

0.725 

0.725 

0.803 

i 

i 

6 

0.645 

0.518 

0.518 

0.784 

1 

7 

0.614 

0.482 

0.482 

0.765 

2 

9 

0.585 

0.447 

0.895 

0.708 

2 

11 

0.501 

0.354 

0.709 

0.689 

2 

13 

0.474 

0.326 

0.653 

0.669 

6 

19 

0.448 

0.300 

1.800 

0.631 

4 

23 

0.399 

0.252 

1  .007 

0.612 

7 

30 

0.375 

0.230 

1.607 

0.593 

1 

31 

0.352 

0.209 

0.209 

0.574 

2 

33 

0.330 

0.189 

0.378 

0.555 

3 

36 

0.308 

0.171 

0.513 

0.536 

1 

37 

0.287 

0.154 

0.154 

0.517 

3 

40 

0.267 

0.138 

0.414 

0.498 

4 

44 

0.248 

0.123 

0.493 

0.479 

2 

46 

0.229 

0.110 

0.219 

0.421 

2 

48 

0.178 

0.0748 

0.150 

0.383 

10 

58 

0.147 

0.0563 

0.563 

0.345 

1 

59 

0.119 

0.0411 

0.0411 

0.326 

6 

65 

0.106 

0.0346 

0.208 

0.307 

8 

73 

0.0941 

0.0289 

0.231 

0.288 

11 

84 

0.0827 

0.0238 

0.262 

0.269 

5 

89 

0.0721 

0.0194 

0.0968 

0.249 

5 

94 

0.0622 

0.0155 

0.0776 

0.230 

10 

104 

0.0530 

0.0122 

0.122 

0.211 

11 

115 

0.0446 

0.00941 

0.104 
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Drop  Diameter, 
D  (mm) 

No.  of 
Drops  n 

Accumulated 
Number,  N 

D2 

(mm)2 

D3 

(mm)3 

nD3 

(mm)3 

0.192 

10 

125 

0.0369 

0.00708 

0.0708 

0.173 

12 

137 

0.0299 

0.00516 

0.0620 

0.154 

11 

148 

0.0236 

0.00363 

0.0399 

0.135 

16 

164 

0.0181 

0.00244 

0.0390 

0.115 

18 

182 

0.0133 

0.00154 

0.0276 

0.0962 

38 

220 

0.00925 

8.90  x  10"4 

0.0338 

0.0770 

37 

257 

0.00593 

4.57  x  10'4 

0.0169 

0.0578 

27 

284 

0.00334 

1.93  x  10"4 

0.00522 

0.0386 

18 

302 

0.00149 

5.74  x  10“5 

0.00103 
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TABLE  2:  COMPARISON  OF  AVERAGED  DIAMETERS 


Technique 

Mass  Mean 
Diam  (mm) 

Mass  Median 

Diam  (mm) 

Regression  Coef  t 
or  Std  Dev  (±mm)+t 

Standard 

0.385 

0.663 

±  1.255  ++ 

Log  N  vs  D 

0.386 

0.781 

0.990  + 

Log  n  vs  D2 

0.534 

0.718 

0.971  t 

UNCLASSIFIED 


CUMULATIVE  NUMBER 


CUMULATIVE  NUMBER 


APPENDIX  A 

Computer  Program  for  Averaging  Drop  Diameters 


non  non 


45  AU.M+1)«A(1.M) 

A ( 2 »V+1 ) *A I  2  »M ) 

IF<m,FO,I)  GO  TO  47 
M*M— 1 
GO  TO  45 
47  A<1.I>*TEMP1 
A(2.I)«TPMP2 
GO  TO  25 

40  A(2#I)*A(2»I)+A(2»J)  ...* 

A ( 2 ♦ J ) *0  *0 
GO  TO  25 
30  CONTINUE 
25  CONTINUE 
IK«0 

00  50  J  * 1 #  N 
IK=IK+1 

IF(A<2.J).FQ.0.0>  60  TO  60 
A!1,IK)«A<1»J) 

A(2«IK)*4(2»J) 

GO  TO  50 
60  IK-IK-1 
5  0  CONTINIJF 
N=- 1 K 

EXPANDING  THE  MATPIX 

SUM=0.0 
SUMl-O.O 
SUM2»0.0 
DO  70  J*  1  »N 

A(3»J)«0.1P2*(A< 1 « J) >**0.997 
A ( 6  » J ) *A ( 3  » J ) **3 
SUM«SUM+A (  2 1 J  > 

A!4,J)»A( 2*J>*A(3» J>**3 
SUM1*SUM1+A(4»J) 

A( 5  » J  > *SUM1 

70  SUM2«SUM2+A<2»J)*A(6.J)**2 
WP I TE ( 6  *  4 )  N 

WRITE! 6,5 >  (All  *  J ) *  A ( 2 • J ) ♦ A ( 3  » J ) ♦ A ( 6  *  J ) *  A ! 4»  J ) . J*1 »N ) 

CALCULATE  MASS  MFAN  AND  ST.  DEV. 

SMFAN*  < SUM1/SUM ) #*0.33333 
AVM AS* SUM l/?,0 
AVMASS*3. 141 56/6. 0*AVMAS 

STDEV* ( SORT ! SUM2/5UM- ( SUM 1/ SUM  >  **2  > >**0.333333 

SUM3*0.0 

00  100  J*1.N 


//AvDROP  JOR  (  ),  'GORDON  HILLS 

//  MSGLFVFL*1 #mSGCLASS  =  A,CLASS-K*TIME*  <0*29 ) 

//•MAIN  LINES* 10 .CARDS- 20 »ORG-RM028 
//•FORMAT  PU»DONAMF=»DFST “LOCAL 
//  EXEC  FORTXCG 
//SYSPRINT  DO  SYSOUT-A 
//SYSTEM  no  SYSOlJT  =  A 
//FORT.SYSIN  DP  * 

C 

C  PROGRAM  TO  EVALUATE  vvq'S 

C 

DIMENSION  A(6«1001> »CIRCLE( 10) 

1  FORMAT ( 2F10 • 3 ) 

2  FORMATCl  INPUT  DATA  FOR  POSITIONS*  SlOAlt'  WHICH  HAS  '♦ 
1 1 4 » •  INPUT  CARDS' ,//,8X»'STAIN  S I Z E • * 5X * « NO.  OF  DROPS**//) 

3  FORMAT (7X»F10.3»5X*F10.1) 

4  FORMAT ( *  1  RFDUCED  DATA  HAS  **I4,*  SETS  OF  DATA'»//.8X. 
l'STAIN  SIZE ' »5X * 'NO.  OF  DROPS ' »5X ,' DROP  SIZE  * »8X* *D**3 * » 
213X*  *NO*»3 • ♦//) 

5  FORMAT (7X,F10«3,5X,F10*1»5X»F10.5»5X»E13*6»5X»E13.6) 

6  FORMAT! 0.'  SOMETHING  WRONG  THE  AVERAGE  MASS-  ••FlO.St'  IS  B 
1THAN  THE  TOTAL  MASS*'  .F10.5) 

7  FORMAT ( ' 1 • »SX  *  *  HERE  ARE  THE  AVERAGED  DATA'.//) 

8  FORMAT ( 1 OA 1 ) 

9  FORMAT ( 5X  ,  •  MEAN- ', F 10 . 4, / »6X STANDARD  DEVI  AT  I  ON- ' *F1C . 6 , / , 
16X  » ' AVFRAGF  MASS*  ' »F10 .4 , / *6X* • MASS  MEDIAN- ' ,F10. 4) 

DO  R98  NJK-1,2 

RFAD (  5  » B  )  !CIRCLF( I )  ,1-1.10) 

DO  10  J-1.1001 
N*J-1 

RFAD (5,1 ) A  I 1,J) ,  A  ( 2  ,  J  ) 

I F ( A(ltJ) .EQ.0.0.AND.A(2»J) *EQ.O.O)GO  TO  20 
10  CONTINUF 

20  WRITE!*,?)  (C1RCLF! I ) ,1-1,10) »N 
WR I TE ( 6 » 3 )  ( A ( 1 ,  J ) « A ( 2  ,  J ) » J* 1 »N ) 

C 

C  DATA  REDUCTION 

C 

DO  75  J-2.N 
IK-J-1 

TEMP1-AH  *J) 

TFMP2-A ( 2  » J ) 

DO  30  I  *  1 » I K 

ROU1 -A (1*I)+0»01*A(1,I > 

ROU2-A(l,I )-O,01*A( 1,1 ) 

I F  <  A ( l.JI.GT.ROUl )  GO  TO  30 

IF(A(1,J).GT, ROU2.AND. A ( 1  * J ) » LT .BOU1 )  GO  TO  40 
M-lK 


BIGGER 


TO  40 


Ann 


100  SUM3»SUM3+.(  (A  (6tJ)~SUVl/SUM2  >**2)*A(2,J) 
ST0EV= ( SORT < SUMS /SUM) >*#0.33333 


CALCULATE  MASS  MEDIAN 


00  «0  J-l.N 

I F( AVMAS.LT.A (5 ♦ J ) >  GO  TO  90 
BO  CONTINUE 

DRMASS«A( 5»N)*3.14156/6.0 
WRITE  ( (it  6 )  A  VMASS  tORMASS 


GO  TO  999 

90  $MED-(AVMAS-A(5»J-1 ) > / ( A < 5 * J > -A ( 5 » J-l ) )* 
1<A(*#J>-A«  3  «  J-l ) )+A(3»J-l> 

WRITE<6»7) 

WR1 TE ( 6  »9 )  $M£AN  »STD£V  .AVMASS. $MED 

998  CONTINUE 

999  STOP 
END 

//GO.SY5LOUT  DO  SYSOUT*A 
//GO.FT06F001  DD  SYSOUT-A 
//GO.FT07P001  DD  SYSOUT-R 
//GO.SYSIN  DO  * 

/* 


} 


//JMDROP  JOA  i  )  *  •  GORDON  HILLS 

//  M$GLFVEL*1 »«SGCLA$S=A»CLASS*K»TIME* (0*29) 

//•MAIN  LINFS*10#CARDS«20»ORG*RM028 
//•FORMAT  PU«nnNAViF  =  f DEST  =  LOCAL 
//  FXEC  FORTXCG 
//SYSPRINT  OD  SYSOUT«A 
//SYSTEM  OD  SYSOUT=A 
//FORT.SYSin  DD  • 

c 

C  PROGRAM  TO  EVALUATE  VMO ' S 

c 

DIMENSION  A(i0*I001)*CIRCLEU2) 

1  FORMAT ( ?F 1 0 • ^ 1 

2  FORMATCl  INPUT  DATA  FOR  POSITIONS*  **12A1**  WHICH  HAS 
114*'  INPUT  CARDS' *//.8X»*STAIN  S I ZE • • 5X » • NO.  OF  DROPS ' ♦ // I 

3  FORMAT (7X*F10.3«SX»F10.1) 

4  FORMAT('l  RFP'  'CED  DATA  HAS  '*I4,»  SETS  OF  DAT  A  '  » /  /  »8X  » 
l'STAIN  SIZE  * *5X» 'NO.  OF  DROPS  * ♦ SX . ' DROP  S I  ZE ' ♦ 10X *  * N ' » 10X * 
2'LOGN* *13X**D*«2' ,//) 

5  FORMAT \jy »F10.3.5X*F10.1.5X»F10.S#5X»F10.1»5X»E13.6»5X»F10.7) 

7  FORMAT ( * i  COEFFICIENTS  FOR  LOGN  VS  D'.//> 

fl  FORMAT ( 1  * A1 ) 

9  FORMAT <  iu  1  INTERCEPT*  ' »F10.6 ♦/ .6X SLOPE*  *»FlC.6»/» 

16X. CRESS  ION  COEFFICIENT*  » .F10.6 »/ »6X ♦ ' M,  MEAN  D.*  '*F10.6* 
2/»6X»'M.  MEDIAN  D.*  ' •F10.6*/// ) 

11  FORMAT  (/'/*•  COEFFICIENTS  FOR  LOGN  VS  D**2'»//) 

DO  Q98  NUX*1»1 

RFAD (5.8)  (CIRCLES  )  .1*1*12) 

no  io  j-i.imi 

N*J-1 

RFAD ( 5 ♦ 1 ) A ( ] »J) «  A  (  2  *  J ) 

IF(A(1 »J) .FO.O.O.AND.A( 2»J) .EO.O.OIGO  TO  20 
10  CONTINUF 

20  WRITER. 2)  (C1RCLFU  )  .1*1*12)  .N 
WP I  TP  <  6  *  3 )  ( A ( 1 ♦ J ) • A ( 2  * J ) *J*1 »N> 

C 

C  DATA  REDUCTION 

C 

DO  25  J*2»N 

TFMP 1 *A ( 1 tJ) 

TFMP2*A ( 2 ♦ J ) 

DO  30  I *1  *  I K 

P0U1*A (1*1 ) +0.01* A (1*1 ) 

ROU2*A(l*I ) “0.01* A l 1*1) 

IF(A(1*J).LT. ROU? )  GO  TO  30 

IF(A( 1»J).0T. BOU? .AND. A ( 1 » J ) • LT .BOUJ )  GO  TO  40 
M*  I X 


u  u  u 


45  A(1»M+1)«A(1»MJ 
A(2»M+1 )«A{ 2»M) 

IFfM.FO.H  <50  TO  47 
MxM-1 

GO  TO  45 
47  A ( 1  *  I ) «TEMP1 
A(2»I»«TFMP2 
GO  TO  25 

40  A(2.I  >»A(2.n^A<2.J) 
A(2.J)»0.0 
GO  TO  25 
30  CONTI NllF 
25  CONTINUE 
I  K*0 

DO  50  J * 1 » N 
IK-IK+l 

IF(A(2»J)»eQ«0«0>  GO  TO  60 
A ( 1 » I K ) *A ( 1 # J ) 

A ( 2  *  I K) »A ( 2  » J ) 

GO  TO  50 
60  I K» I K— 1 
50  CONTINUF 
N«Itf 

FXPANDING  THE  MATRIX 

SUM-0. 0 

SUM1»0.0 

SUM2-0.0 

SUM3«0.0 

SUM4-0.0 

SUM5-0.0 

SUM6«0.0 

SUM7*0.0 

SUM8-0.0 

DO  70  J-l.N 

A(3*J)«0.192*(A( 1  * J ) ) **0.997 
SUM.SUM+AI2.J) 

A(4*J)«SUM 

A ( 5  .  J ) -AL0G19 ( A ( 4 . J ) ) 

A (6  »J ) »A ( 3  » J ) **2 

A(7*J)»A(6»J)**2 

A(fl»J)«A(6»J)*A{5.J) 

A(9»J)«A(3.J)#A<5»J) 

A(10*J)«A(5*J)**2 

SUM1»SUM1+A(6«J) 

SUM2»SUM2+A(7,J) 

SUM3*SUM3+A  <  3  * J  > 


SUM4*SUM4+A(9»J) 

SUM5»SUM5+A(P*J) 

SUM6»SUM6+A(3»J) 

SUM 7* SUM 7+ A (  5  ♦  J  ) 

SUMfl-SUM«+A(10»J) 

70  CONTINUE 

WRI TE ( 6 ,4 )  N 

WRIT£(6»5)  ( A ( I ,  J ) *  A ( 2  *  J ) *  A ( 3  *  J ) » A ( 4, J  ) *  A  (  5  »  J )  •  A ( 6* J ) » J-l »N > 
IF(N.EQ.l)  GO  TO  100 
C 

C  CALCULATE  COEFFICIENTS  FOR  LOG  N  VS  0 

C 

Al* (SU«1*SUM7-SUM3*SUV4) / <N*SUM1-SUM3**2 ) 

Rl» (N*SUM4-SUM3*SUM7 ) / <  N*SUM1-SUM3#*2 > 

Rl« ( N*SUM4-SUM3*SUM7 ) / ( SORT  I  I N*SUM1-SUM3**2 ) * ( N*SUMB-SUM7**2 >  > ) 
DM1-0.789/APSIP1 ) 
n01«l.SR4/ARS(Pl ) 

C 

C  CALCULATE  COEFFICIENTS  FOP  LOG  N  VS  D**2 

C 

A2« ( SUM?*SUM7-SUM1#SUM5 ) / ( N*SUM2-$UM1*#2 ) 
R2«(N*SUM5-SUM1*SUM7)/(N*SUM2-SUM1**2 ) 

R2*  ( N*SUM5-5UM1*5UM7  )  /  (  SORT  (  ( N*SUM2-SUM1  **2  )  *  I  N*SUM8-SUV,7**2  )  )  ) 
DM2-0.724S/SQRT( ABSIP2  )  ) 

D02«0.o733/S0RT( ABS(P2) ) 

GO  TO  150 
100  Al-0.0 
A2»0.0 
Rl-0.0 
R2-0.0 
Rl-0.0 
R2-0.0 
DM1*A(3»1 > 

DM2«A 13*1) 

001*A<3»1> 

D02"A ( 3*1) 

C 

C  PRINT  RESULTS 

c 

150  WR I TF ( A  #  7 ) 

WRITE! A#o)Al#Pl tRlfOMl ,001 
WRITE(ft.ll) 

WRITE«ft#9)A2»P2,R2»DM2,D02 

998  CONTINUE 

999  STOP 
ENO 

//GO.SYSLOUT  00  SYSOUT=A 
//G0,FT06P001  00  SYSOUT-A 


//GO.FT07F001  OD  SYSOUT-B 
//GO.SYSIN  DO  * 
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